Confining stationary light: Dirac dynamics and Klein tunneling 
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We discuss the properties of ID stationary pulses of light in atomic ensemble with electromagnet i- 
cally induced transparency in the limit of tight spatial confinement. When the size of the wavepacket 
becomes comparable or smaller than the absorption length of the medium, it must be described by 
a two-component vector which obeys the one-dimensional two-component Dirac equation with an 
effective mass m* and effective speed of light c*. Then a fundamental lower limit to the spatial 
width in an external potential arises from Klein tunneling and is given by the effective Compton 
length Ac = h/(m*c*). Since c* and m* can be externally controlled and can be made small it is 
possible to observe effects of the relativistic dispersion for rather low energies or correspondingly on 
macroscopic length scales. 

PACS numbers: 42.50.Gy, 42.50.Ct, 41.20.Jb 



When photons are confined to a volume smaller than 
the wavelength cubed their interaction with atoms is 
dominated by quantum effects. This principle has been 
exploited in cavity quantum electrodynamics, where the 
light is confined by means of low-loss micro-resonators 
The tight confinement results in a strong coupling 
which can be used e.g. to build quantum gates between 
photonic qubits. Yet with a decreasing resonator volume 
it becomes more and more difficult to maintain high Q 
values. However, as shown by Bajcsy et al. 0, y], it 
is possible to create spatially confined quasi-stationary 
pulses of light with very low losses without the need of a 
resonator by means of electromagnetically induced trans- 
parency (EIT) 0, [H with counterpropagating control 
fields. For a weak confinement in the longitudinal direc- 
tion, stationary light is well described by a Schrodinger- 
type equation with complex mass for a normal mode of 
the system, the stationary dark-state polariton 0,0, @l- 
We here show that this is no longer the case for stronger 
spatial confinement, i.e. when the characteristic length of 
the photon wavepackets becomes comparable or smaller 
than the absorption length of the medium. Here a de- 
scription in terms of a two-component vector obeying a 
one-dimensional Dirac equation becomes necessary. The 
two characteristic parameter of this equation, the effec- 
tive mass m* , and the effective speed of light c* depend 
on the strength of the EIT control fields and can be made 
many orders of magnitude smaller than the vacuum speed 
of light c and respectively the mass of the atoms forming 
the EIT medium. As a consequence effects of the rela- 
tivistic dispersion can arise already at rather low energy 
scales. On one hand this leads to a fundamental lower 
limit for the spatial confinement of stationary light in an 
external potential given by the effective Compton length 
— h/(m*c*) which due to the smallness of m*c* can 
become large. On the other hand it opens the possibility 
to study relativistic effects such as Klein tunneling Q and 
Zitterbewegung which regained a lot of interest recent! 
in connection with electronic properties of graphene 



and ultra-cold atoms in light- or rotation-induced gauge 
potentials [111 ] . 

We here consider an ensemble of quantum oscillators 
with a double-A structure of dipole transitions as shown 
in FigUl The ground state \g) and the meta- stable state 
\s) are coupled via Raman transitions through the excited 
states \e±) by two counterpropagating control laser of op- 
posite circular polarization and Rabi-frequencies fi± and 
two counterpropagating probe fields E± again of oppo- 
site circular polarization. Both A schemes are assumed to 
be in two-photon resonance with the ground state tran- 
sition, which guarantees EIT. Furthermore the control 
fields are taken homogeneous, constant in time and of 
equal strength fi + = = O = As shown in 

[7[ the control fields generate a quasi-stationary pulse of 
light. For the present discussion we restrict ourselves to 
a one-dimensional dynamical model. One can show that 
the transverse dynamics occurs on a much longer time 
scale than the longitudinal one, which justifies this sim- 
plification [121 ] . 
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FIG. 1: (color online) stationary light scheme: The interac- 
tion of double A atoms driven by two counter-propagating 
control fields of (equal) Rabi frequency Q and opposite cir- 
cular polarization with two counter-propagating probe fields 
E± of corresponding polarization generate a quasi- stationary 
pattern of the probe fields. 



We introduce normalized field amplitudes that 
vary slowly in space and time , E±(r,i) = 

■ ^£±(r, t) exp{— i(ujt =F kz)} + h.a^j , and continuous 
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atomic- flip operators &^(r, t) = 



^GAV(r) With 

& = being the flip operator of the jth atom. 

The sum is over all AN atoms in a small volume AV 
around position r. Then the dynamical equations read 
in the linear response limit, i.e. for a low probe light 
intensity 
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Here we have introduced the common single-photon de- 
tuning of the upper states from the control and probe 

field transitions A. n is the atom density and g = -| yj |^ 
the common coupling constant of both probe fields with 
p denoting the respective dipole matrix element. 7 is 
the transverse decay rate of the optical dipole transitions 
|e±) — \g) and we have used that in the linear response 

limit <Jgg ~ 1. 

Adiabatically eliminating the atomic variables from 
the equations of motion leads to the shortened wave equa- 
tions for the forward and backward (£-) propagat- 
ing field components 



— ±c— 

dt dz 



£± = 



g 2 n 



g 2 n 

2T 



{£+ + £-) 



tan 2 #( j^£ + 



d 

dt 1 



(4) 



where T = 7 + zA, and tan 2 = g 2 n/Q 2 , with ^ 2 = 
Q\ + • If the characteristic length scale of the probe 
fields is large compared to the absorption length of the 
medium, L a b s = 7c/(g 2 n), it is convenient to introduce 
sum and difference normal modes Es = (£++£_)/ a/2 and 
Ed = (£+ — £-)/ V%- Expressing the equations of motion 
in terms of these normal modes and subsequently adia- 
batically eliminating the fast decaying difference normal 
mode Ed leads to a Schrodinger-type equation with a 
complex effective mass 0, [?J . 
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Here m* denotes the real part of the effective mass 
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As comparative scales we have introduced in the second 
equation the mass m and the recoil velocity v rec = hk/m 
of an atom and the group velocity of EIT v gr = ccos 2 0. 

On the other hand if the characteristic length scale of 
the probe fields becomes comparable to the absorption 



length non-adiabatic couplings between the sum and dif- 
ference mode are relevant and the elimination of the dif- 
ference mode is no longer valid. Instead one has to keep 
both amplitudes which can be collected in a two compo- 
nent vector £_ = (£ + ,£_) T . The equation of motion can 
then be written in the compact form 



^ dt ^ dz , 
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where A = 1 + ^ tan 2 0(1 + a x ), B = ca z , and C = 
(a x — l)g 2 n/(2r), a x and a z being the Pauli matrices. 
Applying the transformation 



£ exp {f3a x } E, 



(8) 



with tan(2/3) = (1 - cos 2 0)/(l + cos 2 0) finally yields for 
large single-photon detuning | A| ^> 7, i.e. for a negligible 
imaginary part of the effective mass: 
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Here we have removed an irrelevant constant term by 
a gauge transformation. Eq.(|9|) has the form of a two- 
component, one-dimensional, massive Dirac equation. 
The effective speed of light c* in eq. ([9]) 

c* = c cos 9 = y/v gr c (10) 

which in EIT media can be varied over a large range and 
can be much smaller than the vacuum speed of light. 
It should be noted that despite the formal equivalence 
of equation ([9]) to the two component Dirac equation the 
fundamental quasi-particle excitations of the light matter 
interaction are bosons [7] and can e.g. undergo Bose 
condensation [13]. Equation ([5|) is of course recovered 
from eq.© in the low energy limit of long wavelength 
excitations. 

The characteristic length scale at which relativistic ef- 
fects become important is the Compton length 
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While for electrons Ac is on the order of pico-meters, 
for stationary light it can become rather large. It can 
exceed the absorption length if the EIT group velocity is 
sufficiently large, i.e. v gr /c > 7/ A. Since typical values 
for the optical depth OD = Lj L a b s of EIT media are in 
the range between a few and a few hundred, Ac can be 
a sizable fraction of the medium length L and thus can 
become macroscopic. 

The free evolution of a quasi-stationary pulse of light 
is quite different in the two limits L ^> Ac, eq.©, and 
L ^ Ac, eq.©, which allows for a simple experimen- 
tal distinction of the two regimes. This is illustrated in 
Figl2]which shows the dynamics of a wavepacket obtained 
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from a numerical solution of the full Maxwell- B loch equa- 
tions, which agrees with the dynamics from the effective 
Schrodinger and Dirac equations. One recognizes in the 
first case the familiar slow dispersive spreading, while in 
the second case two split wavepackets emerge which move 
outward with the effective speed of light c* . 




FIG. 2: top: diffusive expansion of an initial gaussian station- 
ary light pulse with width L p = 2.5Ac- bottom: the same for 
an initial width of L p = 0.05Ac- Results are obtained from 
numerical solutions of the full Maxwell-Bloch equations which 
agree very well with solutions of the Schrodinger equation © 
(top plot) and the Dirac equation ((9]) (bottom plot). The pa- 
rameters are 7/ A = 0.01, ft±/A = 0.2, cos (9 = 0.9975 and 
hence \c/L a bs = 199.5 



It is well known that a relativistic wave equation does 
not permit a confinement of a wavepacket to less than 
the Compton length. E.g. a square-well potential 
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with a — > and Uq — > 00 such that Uqcl = const, has 
lowest energy eigensolutions with energy 
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The corresponding eigensolutions have the form 
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Thus the characteristic confinement length L CO nf reads 
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which is always larger than Ac- One can show in general 
that any eigensolution of any ("electrostatic") confining 
potential has a minimum size of 4p. If solutions with an 
energy exceeding ^m*c* 2 exist they are resonances which 
have a finite width due to Klein tunneling into the neg- 
ative (positive) energy continuum [9]. Depending on the 
form of the confining potential the corresponding decay 
rate can however be small. The effect of Klein tunnel- 
ing is illustrated for confined stationary light in FigOl 
Here an initial gaussian stationary wavepacket of initial 
width L p = 0.05A C is considered in a potential well with 
a = 0.1A C and potential depth U = 1.875ra*c* 2 (top) 
and Uq = 3.125m*c* 2 (bottom). Due to the mismatch of 
the initial wavepacket with the bound-states of the poten- 
tial there is some initial loss. After some time the pulse 
shape remains rather unchanged however in the first case 
while it displays continuing decay in the second due to 
Klein tunneling. The plots are obtained from a solution 
of the full ID Maxwell-Bloch equations with an addi- 
tional potential generated by a finite, space-dependent 
two-photon detuning which show very good agreement 
with the solutions of the corresponding Dirac equation. 

There is another well-known effect of the Dirac dy- 
namics which can be observed in the present system. 
If we consider an initial stationary gaussian wavepacket 
and switch the relative phase between the two counter- 
propagating drive fields instantaneously from to 7r/2, 
the relative motion of the two emerging wavepackets is 
superimposed by an oscillation, known as Zitterbewe- 
gung. After the tt/2 phase flip the initial wavepacket 
reads in k space 
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Then in the large-time limit t ^> fi/(m*c* 2 ) one finds 
for the center of mass of the two wavepackets (setting 
7/ A = 0) 
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which shows the characteristic oscillation with frequency 
2m*c* 2 /ft. Fig HI shows the center of mass of a pair of 
left and right moving wavepackets for an initial gaus- 
sian pulse obtained from a numerical solution of the ID 
Maxwell-Bloch equations as well as the analytic result 
([T7]) obtained from the Dirac equation. The Zitterbewe- 
gung can be observed e.g. by detecting the overlap of 
the left and right moving wavepackets after exiting the 
medium. A typical scale of the amplitude of the Zitter- 
bewegung would be L a b s /10 ~ 0.1cm which should be 
much easier to observe than the values in the range from 
l-100nm predicted for cold atoms [ll|, graphene [14], ions 
[lH and photonic crystals [1^ . 
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FIG. 3: (Color online) evolution of an initial gaussian station- 
ary light pulse with initial width L p = 0.05A C in a square well 
potential (the extension of which is indicated by dashed red 
lines) with a = 0.1 Ac and depth Uo = 1.875m* c* 2 (top), and 
Uo — 3.125m* c* 2 (bottom) respectively. Besides initial losses 
due to mismatch of the wavefunctions (see insets at the top) 
a fast decay due to Klein tunneling is apparent in the second 
case. 



In summary we have shown that stationary light in 
the limit of tight longitudinal spatial confinement must 
be described by a two- component, one-dimensional Dirac 
equation with effective mass and effective speed of light 
that can be controlled externally and that can be much 
smaller than the corresponding values for atoms and light 
in vacuum. As a consequence relativistic effects related 
to the Dirac dispersion can be observed at rather low en- 
ergy scales or respectively at rather large length scales. 
One immediate consequence of the latter is the impossi- 
bility to spatially compress a stationary light pulse below 
the Compton length. Moreover in contrast e.g. to elec- 
trons in graphene, interactions between stationary-light 
polaritons are very week. Thus stationary light may be 
employed to observe relativistic phenomena related to 
the Dirac dynamics in the absence of interactions under 
experimentally realistic conditions. 
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FIG. 4: (Color online) left: center of mass of emerging 
left and right moving wavepackets obtaind from solution of 
Maxwell-Bloch equations of initial gaussian wavepacket of 
width Az = lOLabs < Ac = 199.98L a bs (solid blue line) and 
from solution (fT7)) of the Dirac equation (dashed red line). 
At t = a flip of the relative phase of the drive fields was 
assumed corresponding to £+(0) = i£-(0). Parameters are 
as in fig|5] right: Fourier-spectrum of (z(t)) showing the pro- 
nounced peak at 2m*c* 2 //i. 
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